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Abstract 

We calculate two-loop matching conditions for all the operators that are relevant to 
B -^ Xsl^l~ decay in the Standard Model. In effect, we are able to remove the ±16% uncer- 
tainty in the decay spectrum, which was mainly due to the renormalization-scale dependence 
of the top-quark mass. We find 1.46 x 10^^ for the branching ratio integrated in the domain 
m^+i-/ml G [0.05,0.25], for I = e or fi. There remains around 13% perturbative uncertainty 
in this quantity, while the non-perturbative effects are expected to be smaller. 



1 Introduction 

The forthcoming measurement of the inclusive decay mode B — > Xsl^l^ is expected to 
provide an important test of possible new physics effects at the electroweak scale. However, 
the existing theoretical predictions for the branching ratio in the Standard Model (SM) still 
suffer from many uncertainties, some of which are larger than the expected experimental 
errors. 

The most important theoretical uncertainties are due to intermediate cc states. Because 
of the non-perturbative nature of these states, the differential decay spectrum can be only 
roughly estimated when the invariant mass of the lepton pair mf+i- is not significantly below 
rrtji^. It remains questionable whether integrating the decay rate over this domain can 
reduce the theoretical uncertainty below ±20% [|I|. 

On the contrary, for low s = mf+i-fmlp^i^ (accessible to / = e or /x), a relatively precise 
determination of the decay spectrum is possible using perturbative methods only, up to 
calculable HQET corrections. The dominant HQET corrections were evaluated in refs. ^- 
and found to be small (smaller than 6% for 0.05 < s < 0.25). Effects of similar size are 
found in this region when purely perturbative expressions for cc contributions are compared 
with the ones obtained via dispersion relations in the factorization approximation (see fig. |I] 
in section ^). Thus, the B -^ Xgl^l" decay rate integrated over this region of s should be 
perturbatively predictable as precisely as the B — > Xg^y decay rate, i.e. up to about 10% 
uncertainty. 

Unfortunately, the presently available perturbative calculations |^, |^ have not yet reached 
this precision, even though they are performed at the next-to-leading (NLO) order in QCD. 
The formally leading-order term is (quite accidentally) suppressed, which makes it as small 
as some of the NLO contributions. Consequently, some of the formally next-to-next-to- 
leading (NNLO) terms can have an effect larger than 10% on the differential decay rate. It 
can be easily verified by varying the renormalization scale at which the top quark mass is 
renormalized in the formulae of refs. ^ § . 

The formalism of effective theories, which is conventionally used in the analyses of weak B 
decays, allows the identification of three types of NNLO contributions to 5 ^ X^l^l^ . The 
first type originates from two-loop matching between the Standard Model and the effective 



theory amplitudes, i.e. to two-loop contributions to the Wilson coefficients in the effective 
theory at the scale /io ~ My/- The second type is due to the three- loop renormalization group 
evolution of the Wilson coefficients down to the scale /i^ ~ 771^. The third type originates 
from two-loop matrix elements of the effective theory operators between the physical states 
of interest. One should include one- loop Bremsstrahlung corrections as well. Performing a 
complete NNLO calculation is thus a very involved task. 

In the present paper, we shall calculate only the ffist type of corrections, i.e. those 
originating from the two-loop matching conditions. Our results will allow us to remove the 
significant uncertainty of the former NLO prediction stemming from the dependence on the 
scale /io- The remaining uncalculated NNLO effects will be estimated in section |^. 

Our paper is organized as follows. In section 0, we introduce the effective theory and 
present a complete set of the matching conditions up to two loops. The resulting formulae 
for the so-called effective coefficients are given in section ^ Section ^ is devoted to discussing 
phenomenological implications of our results for B — > XJ^l^. Technical details of the 
matching computation are relegated to section ^ There, one can find an extensive description 
of the two-loop matching procedure for the photonic penguin diagrams, which has been the 
most involved original part of our calculation. Section ^ can serve as a practical guide for 
performing any two-loop matching computation, not necessarily in the domain of flavour 
physics. 

2 Summary of the two-loop matching conditions 

The effective theory lagrangian relevant to B —* Xgl^l" decay has the following form 

*--e// = *-'QCDxQed['U, CI, s, c, b, e, ix,t) 

+ ^ KsVubiC'^Pl" + C^2P2) + K:V,t{C',P^ + C^2P2)] 

+ 



-7=-Y.[iKsVut+v:,v,,)c: + v,:v,,ci]p,. 

V^ i=3 



i=3 

For further convenience, we refrain from using unitarity of the CKM matrix V in all the 
analytical formulae here. The ffist term in eq. (|T]) consists of kinetic terms of the light SM 
particles as well as their QCD and QED interactions. The remaining two terms consist of 



AB = —AS = 1 local operators of dimension < 6, built out of those light fields:[] 

P^= {-sl1^J''ul){ul1^T%l), 

P{= {sl1,.T-cl){cl7^T%l), 
P3= {sLl^bL)E,{qrq), 

P5= {sL7^.,j^,i^,bL)E,{qr'r'r'q), ^^^ 

P7= ^m,{sLa>^''bR)F^,, 
P8= \m,{sLa^^''T%ii)Gl,, 
P9= ^{sLl,bL)Eiiiri), 

PlO= ${sLl,bL)T.l{h^l^l). 

where sums over q and / denote sums over all the light quarks and all the leptons, respectively. 
The Wilson coefficients can be perturbatively expanded as follows 

C? = Cf °) + ^Cf ^) + ^Cf (^) + 0{g% Q = c or t. (3) 

Their values are found in the matching procedure, which amounts to requiring equality of 
b —>■ s+ (light particles) Green functions calculated in the effective theory and in the full 
Standard Model, up to ^[(external momenta and light masses)^/M^]. Contributions of 
order g"^"- to each Wilson coefficient originate from ra-loop SM diagrams, which follows from 
the particular convention for powers of gauge couplings in the normalization of our operators. 
Dimensional regularization with fully anticommuting 75 has been used in our matching 
computation. Using this simple scheme could not cause any difficulties, because the choice 
of the four-quark operator basis in eq. (Q) allowed us to avoid the appearance of Dirac traces 
containing 75 in the effective theory diagrams |^. No such traces were present in the SM 
diagrams, either. 



The MS scheme with scale /io ~ M\y was used for all the QCD counterterms, both in 



^ The s-quark mass is neglected here, i.e. it is assumed to be neghgibly small when compared to nii,. Of 
course, no such assumption is made concerning rric or rur- 



the SM and in the effective theory.0 In addition, several non-physical operators had to be 
included on the effective theory side, because the calculation was performed off-shell (see 
section |] and the appendix for details). 

The 't Hooft-Feynman version of the background field gauge was used for all the gauge 
bosons. It allowed us to perform the matching without making use of the CKM-matrix 
unitarity. 

The only relevant off-shell electroweak counterterm (on the SM side) proportional to s0b 
was taken in the MOM scheme, at g^ = for the s0b term, and at vanishing external 
momenta for the terms containing gauge bosons. 

The obtained matching conditions are the following. At the tree level, all the Q vanish, 
except for 

Cf^ = -1. (4) 

The one- and two-loop matching conditions are summarized below: 
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^ The only exceptions were the top-quark-loop contributions to the renormalization of the light-quark and 
gluon wave functions on the SM side. The corresponding terms in the propagators were subtracted in the 
MOM scheme at q^ — 0. In consequence, no top-quark loop contribution remained in the (W-boson)-(light 
quark) effective vertex after renormalization. 
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The integral representations for the functions Li2 and CI2 are as follows: 
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Our matching results for all the Cj^ are new, except for k = 7,8 and 10. In the cases 



k = 7 and fc = 8, we agree with the previously published results |T0|- The k = 10 case has 
already been discussed by us in ref . |]1T[ , and the original calculation |T^ has been corrected 
in ref. |]13| . 



3 The effective coefficients 



Once the matching conditions are found, the Wilson coefficients should be evolved from 
/io ~ Mw to fib ~ rrib, according to the Renormalization Group Equation (RGE) 
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which has the following general solution 
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In the intermediate step of the above equation, Tg denotes ordering of the couphng constants 
such that they increase from right to left. 

The anomalous dimension matrices 7*^ have the following perturbative expansion 
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The one- and two-loop anomalous dimension matrices have already been evaluated in refs. [0, 
§. However, transforming them to the "new" operator basis (|^) is quite non-trivial (see 
ref. 10 for the 6x6 submatrix). In the "new" basis (and in the MS scheme with the 
evanescent operators specified in the appendix), the matrices 'y'^'-"-* and 7'^*^^-' readQ 
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^ Note that the matrices given here correspond to the normaUzation of operators P7, ..., Piq as in eq. (0) 
and to their ordinary Wilson cocfRcients, not to the so-called "effective" ones that will be introduced below. 



where /3o = y and /?i = ^. The analogous matrices 7*'^°-' and 7*'^^-' can be obtained from 
the ones above by removing the first two rows and the first two columns. 

The complete NNLO prediction for BR[B — > Xgl^l^] depends on two entries of ?7'^*^^^(/ife, /io) 
i.e. on U^2 (fJ'b, /^o) and f/92 {fib, /^o) that are generated by the three-loop matrix 7^^*^^). Unfor- 
tunately, the only entries of 7'^*^^-' that have been calculated so far are the ones corresponding 
to the mixing {Pi, ..., Pe} -^ {Pj, Ps} 13 • Therefore, t/72 (/^fc^ /^o) is known but t/92 (/if,, /io) 
is not. Below, we shall include the unknown Ug2 (l^b, f^o) in our analytical formulae. Its po- 
tential numerical relevance will be tested in the next section. 

After performing the RGE evolution, one evaluates the perturbative expression for 
dr[b -^ Xsl^l^]/ds. It amounts to calculating perturbative matrix elements of the oper- 
ators Pi among the external partonic on-shell states, multiplying them by the appropriate 
Wilson coefficients and performing the phase-space integrals. At NLO, one obtains [^ §: 



dT{b -^ Xsl+l-) 



GlmlpoielKVtbl^ Id, 



ds 487r3 



An 



+ 4 1 + 



s?x 



{C!;"f + 12C:;"Re{Cl^^{s)) 



{21 



+ ^T^(Ct'^^ + Sk9ACf" 



The quantities C^-'-^ can be split into top- and light-quark contributions: 
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Calculating the differential decay rate with the help of eq. (]28|), one must retain only terms 
linear in uj{s) and also set u;(s) to zero in the interference term proportional to Re(C9 (s)). 

The coefficients multiplying Cf^,..., Cg in eqs. ( |30| ) and ( [3T| ) are different from the corre- 
sponding ones in refs. [0, §], because we use a different operator basis here. 

Substituting the evolved Wilson coefficients to eqs. (|30D-(p3D, we obtain the following 
expressions for the "effective coefficients" : 
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where r/ = a,(/io)/as(/ib) and a^ = (^, §, ^, -i|, 0.4086, -0.4230, -0.8994, 0.1456, -1) 
The "magic numbers" entering the above expressions are collected in tables 1, 2 and 3. 
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It is straightforward to verify that our results for the 0{l/as) and 0{1) parts of Cg 
and CIq are identical to the ones found in refs. 0, |]. Only the 0{as) parts are new here. 
As far as C^ is concerned, we just reproduce the result of ref. |1^, where the 0{as) part 
was already present. 

In order to obtain the complete NLO prediction for the B — ^ Xsl^l~ decay rate, one should 



i 


1 


2 


3 


4 


5 


6 


7 


8 


hi 

h'1 

hf 

e'\ 


42678 
30253 

4246707584 

400095925 

3344583818789933 

360615755431797 

199164 

30253 

4298158 

816831 


86697 
103460 

89606166 
13682585 
90790555261878016 


_3 

7 

45043984 
9898119 
6473 
7406 
_6 
7 




1 

14 

34505657 

45891279 

9371 

22218 
2 
7 




-0.6494 
2.0040 
-2.7231 
-2.0343 
-1.9043 


-0.0380 
0.7476 
0.4083 
0.1232 

-0.1008 


-0.0186 

-0.5385 

0.1465 

0.1279 

0.1216 


-0.0057 
0.0914 
0.0205 

-0.0064 
0.0183 


13088650734603675 

115596 

25865 

8516 

2217 



Table 1. "Magic numbers" entering the expressions for Cj and C-j . Three- loop 
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anomalous dimensions from ref. [141 have been used in their evaluation. 
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i 


3 


4 


5 


6 


^7 

i 


8 


9 


pf^' 


80 
203 


8 
33 


0.0433 


0.1384 


0.1648 


-0.0073 


4704688 
25088393 


rt 


3085 
3703 


129 
1058 


-0.1642 


0.0793 


-0.0451 


-0.1638 





' i 


64730 
322161 


18742 
52371 


0.0454 


-0.3719 


-0.3254 


0.0066 


1775737 
809303 


' i 


40 
203 


8 
33 


0.0339 


-0.1122 


-0.2841 


-0.0020 


27051956 
75265179 


' i 


20 
609 


4 
99 


-0.0021 


0.0289 


0.0174 


0.0010 


8908520 
75265179 




316900 
299943 


51388 
128547 


1.9957 


-0.8153 


0.1488 


-0.2353 





< 


183859 
42849 


130739 
128547 


-0.0939 


-0.9763 


0.0393 


-2.2799 





' i 


8129495 
5798898 


4447705 
942678 


0.6261 


-3.6869 


0.2246 


0.0121 


4896690443 
677386611 


r't 


6170 
3703 


258 
529 


-0.5145 


-0.2571 


0.3111 


-0.1829 





' i 


97850 
33327 


398258 
157113 


0.6618 


-2.2108 


-1.6839 


0.0472 


4704688 
2595351 


' i 


20 
21 


4 
9 


0.1833 


-0.2481 


-0.1096 


-0.0090 





' i 


20 
63 


4 

27 


0.0611 


-0.0827 


-0.0365 


-0.0030 





' i 


87527 
99981 


6217 
85698 


-0.1685 


0.0323 


-0.0475 


-0.2018 





Si 


40 
21 


4 
9 


0.2340 


0.3061 


0.0636 


-0.0322 





s'f-' 


1373012 
128547 


735748 
385641 


2.1605 


0.3356 


0.8434 


-0.2456 





s'- 


129460 
33327 


37484 
14283 


0.9813 


-3.2900 


-0.5020 


0.1151 





s'f 


80 
21 


16 
9 


0.7330 


-0.9925 


-0.4383 


-0.0359 





n 


12 

7 


2 
3 


0.1658 


-0.2407 


-0.0717 


0.0990 





^ i 


33606 
3703 


6046 
4761 


-0.1681 


1.2986 


-0.3397 


0.4766 





n 


12946 
3703 


18742 
4761 


0.6951 


2.5871 


0.5664 


-0.3540 





,icL 

'' i 


24 

7 


8 
3 


0.5193 


0.7805 


0.4945 


0.1106 





ut 


2 

7 





-0.2559 


0.0083 


0.0180 


-0.0562 





u'f-' 


168155 
99981 


166 

81 


-1.0892 


-1.1627 


-0.2197 


-0.2193 





< 


6473 
11109 





-1.0733 


-0.0897 


-0.1424 


0.2008 





icL 


4 

7 





-0.8018 


-0.0271 


-0.1243 


-0.0627 





< 


1 

7 


1 

6 


-0.1731 


-0.1120 


-0.0178 


-0.0067 





W'f-' 


251737 
199962 


117137 
85698 


-1.1732 


-0.5134 


-0.3895 


0.0190 





< 


6473 
22218 


9371 
9522 


-0.7257 


1.2038 


0.1408 


0.0238 





,cL 


2 

7 


2 
3 


-0.5421 


0.3632 


0.1229 


-0.0074 






Table 2. "Magic numbers" entering the expression for C; 
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i 


5 


6 


7 


8 


9 


4'^' 


0.0318 


0.0918 


-0.2700 


0.0059 


33160 
235941 


r'\ 


-0.4817 


0.2104 


0.2956 


0.5246 





' i 


0.2164 


-0.4330 


-0.9126 


0.0660 


6672596 
12976755 


s'\ 


0.6862 


0.8125 


-0.4165 


0.1031 





t'\ 


0.4861 


-0.6389 


0.4699 


-0.3171 





u'\ 


-0.7505 


0.0221 


-0.1182 


0.1799 





< 


-0.5075 


-0.2973 


0.1168 


0.0213 





vT' 


-0.1242 


-0.0956 


-0.1628 


-0.0176 


157366 
393235 



Table 3. "Magic numbers" entering the expression for C^^ . 

use eqs. (|28|)-(|33D and neglect the 0{as) contributions to the effective coefficients Cj:^{s) 
(i.e. include only the 0{l/as) and 0{1) parts of them). On the other hand, in the complete 
NNLO calculation, it is not sufficient to take into account the 0{as) parts of the effective 
coefficients. One should also modify eq. (|28D by including effects originating e.g. from two- 
loop matrix elements of the four-quark operators and the corresponding Bremsstrahlung 
corrections. 

In the present paper, we are able to include the NNLO effects only partly. We shall simply 
use eq. (pSD, but at the same time we will include the 0{as) contributions to the effective 
coefficients. In this way, we will include all the m-f-dependent NNLO contributions to the 
branching ratio,Q as well as the terms enhanced by 1/s^ ~ 4.3. It is important to calculate 
the mt-dependent terms at the NNLO level, because both Cg{fio) and C^o(/^o) grow with rrit 
in the formal limit rrit —>■ oo. Therefore, mf/M^^ ~ 4.8 plays the role of an enhancement 
factor, too. 

Above, we have presented explicitly all the 0{as) parts of the effective coefficients. How- 
ever, the unknown quantity Uq2 (yUfe,yUo) occurred in C'g^-'-' (s). In our numerical calculations 
described in the next section, it will be assumed that f/92 (fit, fj^o) vanishes. We shall relax 
this assumption below eq. (^91) , and check that the expected numerical effect of f/92 (/i^, /io) 
on the decay rate is very small. 



"* The only exceptions are the mt-dependent contributions from the one-loop matrix elements of P7 and 
Pg. However, they are proportional to the relatively small Wilson coefficients Crdib) and Csil^b) that do not 
grow with rrit in the formal limit nit ^ 00. 
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4 Phenomenological implications 

In the present section, we shall study the numerical importance of the calculated NNLO 
effects as well as the uncertainties due to the yet unknown contributions. 

As a first step, let us calculate the effective coefficients for several different values of 
fio and fify. We will vary fif, by a factor of 2 around mi, ~ 5 GeV, i.e. we will take 
fih = 2.5, 5 and 10 GeV. In the expressions for Cl'^ and ACg, we will vary /iq by a 
factor of 2 around Mw ~ 80 GeV, i.e. we will take /io = 40, 80 and 160 GeV. In the expres- 
sions for Cf^^ , we will vary /io by a factor of 2 around y/Mwmt ~ 120 GeV, i.e. we will 
take /io = 60, 120 and 240 GeV. 



The remaining input parameters will be equal to |15 



pole _ 1 VQ Q n^\r i\/r... _ on /ii n^\r ,.2 



as{Mz) = 0.119, mr^ = 173.8 GeV, Mw = 80.41 GeV, s; = 0.23124. 



'e// r^Qeff ^^1 J\n^ff 



Since we shall keep s arbitrary, our expressions for Cg , Ciq and ACg will read 

m 



CQeff ^ A'i + I^u{s)+T^h[-^^,s\+U^h{l,s) + W^h{S),s), (42) 

C%'^^ = A% + R%ujis), (43) 

(44) 



Act" 



'9 



/i(0,s)-/iP|,s 



mt 



''b 

The coefficients A^ , ..., W^ are independent of rric, and they only weakly depend on nib via 
the logarithm ln(r7ife//ife). In this logarithm, we shall use nib = 4.8 GeV. 

In tables 4 and 5, our results for C^ , A^ , ..., W^ and Zg are given, both with and 
without the 0{as) contributions. They allow the following observations: 

• The dominant contributions to the "effective coefficients" and to the decay rate origi- 
nate from Aq and A\q. However, the coefficients Cj ^ are not much less important, 
because of the factor "12" in the last term of eq. (pHf). 



The inclusion of the 0{as) contributions significantly reduces the /io- dependence. It is 
especially important in the case of A\q, which had varied by more than ±10% before 
including the 0{as) correction. The dependence on /io remains significant only in the 
relatively small quantities such as R\q. {R\q is multiplied by uj{s) G [—1.32, —1.24] for 
se [0.05,0.25]). 

13 



/io [GeV] 
fib [GeV] 


40 

5 


80 
5 


160 

5 


80 
2.5 


80 
10 


asifio) 

Oisifib) 

V 


0.136 
0.215 
0.633 


0.121 
0.215 
0.565 


0.110 
0.215 
0.510 


0.121 
0.267 
0.454 


0.121 
0.180 
0.674 


Cr^^ withC(a,) 
^7'^^ without 0{as) 


0.567 
0.631 


0.567 
0.631 


0.566 
0.632 


0.554 
0.634 


0.579 
0.631 


Al with 0{as) 
A% without 0(a,) 
A% only 0(l/a,) 


-4.685 
-4.620 
-1.569 


-4.683 
-4.635 
-1.964 


-4.689 
-4.681 
-2.315 


-4.612 
-4.750 
-2.181 


-4.828 
-4.635 
-1.612 


B^ with 0(a,) 
i?g without C(a,) 


-0.315 
-0.107 


-0.316 
-0.134 


-0.320 

-0.158 


-0.415 
-0.186 


-0.242 
-0.092 


Tg^ with 0(a,) 
Tg^ without C(a,) 


-0.641 
-0.505 


-0.625 
-0.374 


-0.603 
-0.255 


-0.393 
-0.115 


-0.807 
-0.576 


[/g^ with C(a,) 
\Jl without 0(a,) 


-0.048 
-0.026 


-0.050 
-0.032 


-0.052 
-0.038 


-0.070 
-0.045 


-0.035 
-0.022 


VTg^ with 0(a,) 
VTg^ without 0(as) 


-0.045 
-0.026 


-0.046 
-0.032 


-0.047 
-0.038 


-0.062 
-0.044 


-0.033 
-0.022 


A^o with 0[a,) 
A\q without 0(a,) 


1.128 
1.081 


1.123 
1.081 


1.119 
1.081 


1.123 
1.081 


1.123 
1.081 


R\^ with 0(a,) 
i?5o without 0{a,) 


0.074 



0.074 



0.074 



0.092 



0.062 



Zg with C(a,) 
Zg without C(as) 


-0.648 
-0.506 


-0.634 
-0.376 


-0.613 
-0.257 


-0.410 
-0.118 


-0.811 
-0.577 



^ce// 



Table 4. C-j , A\, ..., W^ and Zg for various values of fiQ and fib- 



The dependence on fi^ remains rather strong in most of the listed quantities. It follows 
mainly from the fact that two-loop matrix elements of the four-quark operators have 
not been included. It is relevant especially to the cases of C\^^, Tg and Rg, which will 
cause considerable /Xb-dependence of the final prediction for the decay rate. 

The coefficients Wg turn out to be very small, while ACg in eq. ( ^91) is multiplied by 
\(yus^ub) / {Vts^tb)\ — 0.08. In consequence, the terms containing h{0,s) contribute by 
less than 3% to the differential decay rate for s > 0.05, because \h{0, s)\ = |^ — |(lns — 
^7^)| is smaller than 2.2 in this region. This is fortunate, because h{0,s) is expected 
to receive huge non-perturbative contributions from intermediate light hadron states.^ 
The smallness of Wg and V^b allows us to use only the perturbative expression for 
h{0, s) below. We could equivalently just neglect it. 



^ These contributions are expected to be of the same size as h{0, s) itself, after taking an average over a 
sufficiently wide region of s. 
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/io [GeV] 
^^b [GeV] 


60 

5 


120 

5 


240 
5 


120 

2.5 


120 
10 


mf ^(/io) [GeV] 


180 
0.127 
0.215 
0.591 


170 
0.114 
0.215 
0.531 


162 
0.104 
0.215 
0.483 


170 
0.114 
0.267 
0.427 


170 
0.114 
0.180 
0.635 


^7'^^ withO(a,) 
6-7 '-^^ without 0{as) 


0.261 
0.325 


0.265 
0.310 


0.266 
0.297 


0.225 
0.274 


0.300 
0.344 


A\ with 0{as) 
A\ without 0{as) 


-0.547 
-0.425 


-0.541 
-0.506 


-0.544 
-0.579 


-0.541 
-0.509 


-0.542 
-0.504 


Rl with 0{as) 
Rl without 0{as) 


-0.029 



-0.035 



-0.040 



-0.043 



-0.029 



Tg* with 0{as) 
Tg* without 0{as) 


0.0002 



0.0003 



0.0004 



0.0005 



0.0001 



Ui with 0(a,) 
f/^ without 0{as) 


-0.002 



-0.002 



-0.002 



-0.002 



-0.002 



Wl with 0(a,) 
Wl without C(a,) 


-0.002 



-0.002 



-0.002 



-0.002 



-0.002 



A*io with 0{as) 
A\q without 0{as) 


-3.051 
-3.688 


-3.115 
-3.292 


-3.107 
-2.964 


-3.115 
-3.292 


-3.115 
-3.292 


R{q with 0{as) 
R{q without 0(aj 


-0.252 



-0.225 



-0.203 



-0.280 



-0.189 




yteff 



Table 5. C-j and Ai, ..., Wl for various values of /iq and jjb- 



Huge non-perturbative contributions occur in h{m^^/in1^s) as well, for s > {2'mc/mf,)'^. 
It is illustrated in fig. [l|. Dashed lines show the real and imaginary parts of h{z, s) from 
eq. (0), with z = (1.4/4.8)^ and with h{z, 0) subtracted. Solid lines present non-perturbative 
estimates of the same quantities obtained using the formulae and parameters from ref. [^ 
where the factorization approximation and dispersion relations were used.0 

While the solid lines in fig. [1| should not be regarded as the true non-perturbative results 
(because of the factorization approximation), they give us qualitative information on the 
size of expected non-perturbative effects. In particular, we can observe that replacing the 
solid lines by the dashed ones in the region s G [0.05, 0.25] should have quite a small effect 
on the predicted differential decay rate, owing to the relatively small size of Tg in tables 4 
and 5. Actually, the yU;,- dependence of Tg is numerically more important. Our aim below 
will be predicting the decay rate integrated over s from 0.05 to 0.25. We shall use the purely 

perturbative expression for h{z, s), keeping in mind that the /if,-dependence of our prediction 

^ However 4m|, is replaced by 4m^ in eq. (3.4) of ref. [p6|. We thank F. Kriiger for confirming that this 
was a misprint. 
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/ i 




Figure 1: Perturbative and non-perturbative versions of Iie[h{ml/ml, s) — h{ml/ml, 0)] and 
Im[/i(m^/m^, s) — h{ml/ml, 0)] as functions of s (see the text). 

is expected to be larger than the uncertainty stemming from neglected non-perturbative 
effects/ 

As far as h{l, s) is concerned, the argument for using the purely perturbative expression 
can be the same as for h{0, s) (small coefficients) or the same as for h{m1/m1, s) (convergence 
of the perturbative and non-perturbative results for small s) . 



The decay rate given in eq. (^) suffers from large uncertainties due to rnlp^i^ and the 
CKM angles. One can get rid of them by normalizing to the semileptonic decay rate of the 
6-quark 



2 ™5 



dm 



T[b - X^eu,] = ^^^^^^;^\V^\'9 ( ^ > - ' ""^ 



1927r3 



, m 



b,pole , 



mt 



1 I ' 

ml 



where 



g{z) = I - 8z + Sz^ - z^ - Uz'^lnz 



is the phase-space factor, and 



k{z) 



2as{mb) h{z) 



(45) 



(46) 



(47) 



Svr g{z) 

is a sizeable next-to-leading order QCD correction to the semileptonic decay [T^. The 
function h{z) has been given analytically in ref. [|l^]: 

9N /25 239 25 ,\ , / 4 . 17 o\ 0,0, . 

h{z) = -{1-z^) (x~^^+T^ ) +'2ln2f20 + 902--z2 + — zM +z^ln^z {36 + z^ 

+ (1 - z^) (—-—z + —zA ln(l -z)- 4(1 + 30^2 + z^) In z ln(l - z) 
\ o o o / 



^ The non-perturbative effects estimated in fig. n^ are not included in the HQET correction we shall take 
into account later. 
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-(1 + 16z^ + z'^)[6U2{z) - n^] - 32z^/\l + z) 



n^ - 4:U2i^/z) + AU2{-^/^) - 2 111 Z 111 



1 + A/i. 



Thus, the final perturbative quantity we consider is the ratio 

1 d 



R. 



.1+1- 

quark 



[s] 



r[6 ^ Xf-eUf^] ds 



T{b -^ X,l+l- 



(4J 



'" Rquarki'V [^^ 





075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 



0.075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 



Figure 2: Reduction of /^o-dependence of RiJarki-^)- 



-'Kuarkis) [10-^] 




0.075 0.1 0.125 0.15 0.175 0.2 0.225 0.25 



jl+l- 



Figure 3: Remaining /Xb-dependence of RquarkV^)- 
Our results for RgJarki^) i^ ^^e domain s G [0.05, 0.25] are presented in figs. |^ and ^j. In 

rribVoJE) = lis and \V;,Vtb/Vcb\ = 0.976. The 



Of, 



their evaluation, we have used Ogm 
quantity ACg that is multiplied by Vub has been neglected. The dashed lines represent the 
pure NLO results, i.e. the ones with neglected 0{as) parts of the effective coefficients. The 
solid lines are obtained after including the 0{as) terms. Some of them overlap, and look 
like thick lines. 

In both plots of fig. ||, /i;, = 5 GeV, and three different values of /io are chosen. The left 
plot corresponds to varying /iq by a factor of 2 around ^/M^Pm't in CJ. (as in the first 
three columns of table 5) and keeping it fixed to M]y in C^^ . The right plot corresponds 
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to varying /iq by a factor of 2 around M\y in C^ ^^^ (as in the first three columns of table 4) 
and keeping it fixed to \/Mwrnt in (7^ ^^^ . 

The importance of including the two-loop matching conditions is clearly seen: the de- 
pendence on hq decreases from ±16% to around ±2.5% at the representative point s = 0.2. 
Most of the effect is due to the strong mj-dependence of A\q and to the /iQ-dependence of 
mp(//o). 

In fig. D, the scale /iq is fixed to 120 GeV in Cj^ and to 80 GeV in C^^ , while the 
scale /ift takes the values of 2.5, 5 and 10 GeV. One can see that the /ift-dependence increases 
after taking into account the 0{as) contributions to the effective coefficients. When the 
0{as) terms are not included, an accidental cancellation of the /Xfe- dependence occurs among 
the four contributions to the differential decay rate in eq. (p8D. This cancellation becomes 
exact at s ~ 0.06. The 0{as) term that plays the major role in changing the /if,-dependence 
of Aq (see table 4) and in removing this cancellation is proportional to the product of 
Ci (/io) = —15 — 6L from the matching conditions and Inmh/nb from the one-loop 
matrix element of P{^. A future calculation of the two-loop b -^ sl^l^ matrix elements of the 
four-quark operators is desirable, because it should significantly reduce the /it-dependence 
of the prediction for R^Jarki-^) ■ 

When the results described by the solid lines in fig. ^ are integrated over s, we obtain 

rO.25 

'0.05 

c(2) 

where only the error from /i^- dependence is taken into account. Varying 1/^2 from —10 to 
10 (as promised at the end of the previous section) would increase the uncertainty by only 
0.03. Thus, calculating the three-loop anomalous dimensions in the future is not expected 
to have an important impact on the numerical prediction. 

In the end, we relate the integrand of R^Jarki^) ^o ^^e physically measurable quantity 

rO.25 
/0.05 

0.104[(1.36 ± 0.18) - 0.02 ± 0.06] x 10"^ = (1.46 ± 0.19) x 10^^ (50) 



/■0.25 

/ ds R[Zk{s) = (1.36 ± 0.18) X 10-^ (49) 

J0.05 



(•0.25 I- 

5i?[fi ^ X,/+r],e[o.o5,o.25] = BR[B^X,eu]j ds[R[^;,^{s) + 5y^2R{s) + dy^.R^s 



where, again, only the error from the /it-dependence of R^uarki^) is included. The non- 
perturbative HQET corrections 6i/m'^R{s) and 61/^2 R[s) have been found with the help of 
eq. (32) in ref. p and eq. (18) in ref. P, respectively. The O^l/mf) effects are completely 



18 



negligible for s < 0.25 |T^. The experimental value of 0.104 for the semileptonic branching 



ratio is taken from ref. [15 



It is worth indicating that additional non-perturbative corrections due to the motion of 
the 6-quark inside the B-meson would occur if we wanted to impose additional cuts on the 
emitted lepton energies |2^. Such corrections are absent only when the kinematical cut is 



imposed on nothing but the invariant mass of the lepton pair. 

Of course, translating the restriction s G [0.05, 0.25] to bounds in GeV on the lepton 
invariant mass introduces an additional uncertainty due to the numerical value of mhpoie- 
Since the s-spectrum is almost flat in the considered domain, this additional uncertainty (in 
per cent) will be close to ^crmtpoie/''^b,poie, i-e. rather small. 

Finally let us note that restricting the studied domain of s to [0.05, 0.25] makes the 
integrated B -^ XJ^l^ branching ratio smaller, but at the same time more sensitive to the 
sign of C^ (jjb), when compared to the so-called "non-resonant BR" considered for instance 
in ref. [^. If we changed the sign of C^ •'(/if,), the last result in eq. (^D|) would change 
to 2.92 X 10^^. Thus, extensions of the SM that predict opposite sign of C^^^b) (like 
the MSSM in certain dark-matter-favoured regions of its parameter space) might be tested 
with the help of the integrated BR itself, without considering forward-backward or energy 
asymmetries. 

At this point, we finish our phenomenological discussion, and proceed to describing tech- 
nical details of the two-loop matching computation in the next section. 

5 Two-loop matching for photonic A5 = — AS* = 1 
penguins in the Standard Model 

5.1. Preliminaries 

For processes taking place at energy scales much lower than Myj/, the Standard Model 
can be replaced by an effective theory built out of only light SM fields, i.e. the ones that 
are much lighter that the W-boson. Our goal here is to find two-loop QCD contributions 
to the Wilson coefficients of certain operators in the effective theory. The operators we are 
interested in are the ones giving leading electroweak contributions to the /\B = —AS = 1 
transitions accompanied by either a real photon or a lepton pair emission. In the latter case. 
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we restrict ourselves to processes mediated by a virtual photon, i.e. we do not consider in 
this section the SM diagrams where the W or Z boson couple directly to the lepton line. 

The simplest way to find the Wilson coefficients is to require equality of the off-shell IPI 
amputated Green functions calculated in the full SM and in the effective theory. Up to one 
loop, we need to consider the b -^ 57, 6 — > s gluon and h — > sec functions. At two loops, 
only the b ^ s'y function is necessary. In the cases of 6 ^ 57 and b —^ s gluon, we work 
at the leading order in Oem and up to 0[(external momenta) ^/M^]. In the b — >• sec case, 
external momenta can be neglected. 

We set all the light particle masses to zero in the whole calculation. An exception is the 
6-quark mass, which is being included up to linear order. This means that we maintain rrih 
only in Yukawa couplings and in the 6-quark propagator numerators. The terms of order ml 
are neglected. One can justify this procedure by formally treating the 6-quark mass term as 
an interaction with an external scalar field. 

In addition, all the Feynman integrands are expanded in external momenta before per- 
forming loop integration. Such an expansion, as well as setting all the light masses to zero, 
creates spurious infrared divergences that we regularize dimensionally. As we shall see, all 
these divergences cancel out in the matching conditions relating the full and the effective 
theory Green functions. 

The Feynman integrands for the one- and two-loop Feynman diagrams are generated 
with the help of the program FeynArts [^. After Taylor expansion in external momenta 
and factorizing them out, the integrals remain dependent only on loop momenta and two 
heavy masses: Mw and rrit. Subsequent application of the partial fraction decomposition 

(51) 



(g2 — ?7if)(g2 — ml) m^f — m,2 



q^ — r/if g2 — r/ig 



allows a reduction of all the integrals to those in which a single mass parameter occurs in the 
propagator denominators together with a given loop momentum. Finally, after reduction of 
tensor integrals to scalar ones, the non- vanishing integrals obtained at one and two loops are 
respectively as follows: 

Ci" = _ I ' , , . /t-^-L. (52) 



7r2-^r(l + e) 7 (g2_m2)n^ 

^nin2n3 „4-'?f "P^i i ^\'?. I r„2 ™2\n,, /„2 ™2Nn.,r/„ „ N-Jlr?., ' V^"^/ 
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^4-2er(l + e)2 J (g2_^2)ni(g2_^2)n2[(^^_^2)2]"3' 



with arbitrary integer powers n, ni, n2 and ris, and with m, mi ^ 0. The chosen normal- 
ization makes the results free of trivial common factors. 

In eq. (^3]) we have already made use of the fact that our two-loop scalar integrals always 
have at least one massless term in their denominators. This turns out to be true in all the 
Feynman diagrams we have to consider, provided all the light particle masses are set to zero. 
Therefore, all our two-loop integrals are relatively simple. 

The result for the one-loop scalar integral is 



C(i) 



-D" 



:i + e)n- 



(n-1)! 
which vanishes for n < 0. Here, {a)k denotes the Pochhammer symbol equal to 

a(a + l)(a + 2)...(a + A; - 1), k > I, 
{a)k 



T{a + k) 



r a) 



0, 



(54) 



(55) 



l/[{a-l){a-2)...{a-\k\)], k < -1, 
for integer k and complex a. 

The two-loop integrals can easily be found with the help of Feynman parametrization in 

the cases when mi = m,2 or m,2 = 

2 — ^J-nsi^ + e)ni+n3-3(l + f)n2+n3-3 



^nin2n3 



c(2) = 

nin2?i3 n 

^ '^ O ^^^ _Q 



\ni+n2+n.3+l 



I' 1 \ni+n2+na + 



{ni - l)\{n2 - l)!(ni + ris + ris - 4 + 2e)^, ' 

1 (1 + 2e)n^+n3+n3-5(l + e)n2+n:i~3{^ — ^)l~n2{^ ~ ^)l-na 

(m - l)!(n2 - l)!(n3 - 1)!(1 - e)(l - ivr^eS + 0{e^)) 



(56) 
(57) 



It remains to discuss the case when rrii ^ m,2 and none of the two masses vanishes. The 

(2) 

starting point is the integral ClU, which equals: 



^(2) 



1 



2(i_e)(l-2e) 



1 + x 2 , 

— + -xmx 

e^ e 



;i-2x)ln^x + 2(1 - a;)Li2 ( 1 - -] + 0{e) 

(58) 



where x = m\/m,\ [^. All the integrals with three positive indices can be derived from the 



above result with the help of the following recurrence relations |23|: 



C, 



(2) 
(ni+l)n2n3 



r^(2) 

ni(n2+l)n3 



^nin2(n3+l) 



ni(l- 


-x) 


1 




n2x{l- 


-x) 


1 





[[A-2e-ni-n2-nz + x{ni- n^)]C^^}n 



+ 



xn2 



.^(ni-l)(n2+l)n3 



^ni(n2+l)(n3-l) 



{ [x(4 -2e-ni-n2- n^) + rig - n^]C^}n 



+ 



rii 



.^(ni+l)(n2-l)n3 



'-'(ni+l)n2(n3-l) 



n3(l-x)2 



{[(1 + x)(-4 + 2e) + 2n2 + (1 + 3x)n3]C(2) 
+ 



n2n3 



o r^(2) ^(2) 



713 
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'-'ni(n2-l)(n3+l) 



'-^(ni-l)n2(n3+l) 



} , (59) 



}■ 



All the two-loop integrals defined in eq. ( ^3] ) vanish when either rii or n2 is non-positive. 
When these two indices are positive but n^ is non-positive, they reduce to products of one- 
loop tensor integrals. It is sensible to make this reduction only in the case when the two 
masses are different and non-vanishing. Then we obtain 

r ^3 1 

'~ 2 I -"3-2*; / \ / „, \ ^2-n2+k+j-e 



ci'l 



nmius ^^<(, 



E E 

k=0 j=Q 



-"3 
2fc 



— ?i3— 2fc \ X 
J 



(^_]^^ni+n2+n3+l(^2A;)! 



(ni-l)!(n2-l)!A;!(2-e), 



X 



^(2 ~ e)j+A:(2 — e)_„3_fc^j(l + e)„2_fe_j„3(l + e)ni+na+k+j~3- 



(60) 



Otherwise, one can use eqs. ( ^6]) and (^7]), which apply for non-positive n^, too. Equation 
(|57| ) gives zero in such a case, but eq. ([56| ) does not. 



5.2. The Standard Model side 

Let us start with calculating the b —>■ sj function up to two loops. There is no tree-level 
contribution to this function in the Standard Model. The four IPI diagrams arising at one 
loop are presented in fig. 11. 
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Figure 4: One- loop IPI diagrams for b —>■ s'~f in the SM. The charged would-be Goldstone 
boson is denoted by tt^. There is no W^n^'-f coupling in the background-field gauge. 

We calculate the corresponding unrenormalized amputated Green function off shell, in 
the background-field version of the 't Hooft-Feynman gauge. The Feynman integrands are 
expanded up to the second order in external momenta and mj, (neglecting ml though). As 
in section |^, we refrain from using unitarity of the CKM matrix here. The result can be 
written in the following form: 

•4^^ ^^^.iV«|(K*,K. + K:K.)E/^r5, + V*Vt,f:f^'\x)SA + Oie% (61) 



v^(4 



71 



i=i 



i=i 



where P/j = 1(1 + 75), N^^^ = 1 - e^ + e\^7i'^ + ^k^), k = -f^ -ln{47r) + ln{M^/fil) and 
Sk stand for Dirac structures that depend on the incoming 6-quark momentum p and on the 
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outgoing photon momentum k 



Sj = (tm^^, lfi{p-k), -f^p"^, 7/.^^ pkf,, pPf,, ^p^, j^kf,, 

m6^7^, mb7^^, mbfjf,, nii^'^^p, M'^li?} . ■ (62) 



As we shall see later, explicit results are needed only for the coefficients at the structures 5*2, 
5*8 and Siq. We find 
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X 



(63) 



where a; = rri^/M^. 

Let us now proceed to an evaluation of the first QCD correction to the considered Green 
function. The corresponding two-loop diagrams are shown in fig. |^. 





Figure 5: Two-loop IPI diagrams for b —>■ s'j in the SM. The wavy lines denote either the 
W^-boson or the charged would-be Goldstone boson. The external photon can couple at any 
of the places marked by small circles. 

In analogy to eq. (]6T|), we write the unrenormalized two-loop result as 
^^^A^?M(CK. + K:K.)i:/^f5, + V,:V,,f:fP{x)sA + 0{e), (64) 
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where g is the QCD gauge couphng and A^^^^^ = 1 — 2eK + e^(|7r^ + 2k'^). The two-loop 
analogues of the coefficients given in eq. (|63D are found to have the following form: 
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The last two elements we need to know on the SM side are the b ^ s gluon and h — > sec 
functions up to one loop. They are used to recover one-loop contributions to certain Wilson 
coefficients which take part in the two-loop b ^ s'~f matching condition. 




ti, C, tX \u, C, t M, C, tX \ M, C, t 



b W^ s b TT^ s 

Figure 6: One-loop IPI diagrams for b -^ s gluon in the SM. 




Similarly to the b -^ s'~f case, there is no tree-level contribution to the b ^ s gluon Green 
function in the SM. The one-loop contribution is given by the two diagrams presented in 
fig. p. In analogy to eq. (|6T]), the result can be written as 



^^^5Civ(^) J {y:,y^, + v;:k.) E ufs, + v:^,, i: vf\x)s, ) + o{^), (66) 
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where T" denotes the SU(3) generator corresponding to the outgoing gluon. The coefficients 
at the structures 5*2 , 5*8 and S'lo read 
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Figure 7: Tree-level b —>■ sec diagram on the SM side. 
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Contrary to the functions considered so far, the b — > sec function does acquire a tree-level 
contribution in the SM. It is given by the diagram shown in fig. ^ For vanishing external 
momenta, it give^ 

(69) 



-^^V:,V,t{l,PL)^{YPL)■ 



V2 




Figure 8: One- loop b — » see diagrams on the SM side, which do not vanish in dimensional 
regularization when all the light particle masses are set to zero. 



The non-vanishing one-loop diagrams for the b —>■ sec functions are shown in fig. |^. When 
the external momenta are set to zero, we find the following result for the corresponding 
amputated Green function: 

1 3 



V2 {Any '^'^'"'^ 



6 3Q 

- - 15 - ye ) {-/.PlT-) ® {YPlT^ 



+ 0{e) X 



^ The tensor product symbol F (g) F' is used here to denote the tree- level (sFc)(cF'6) amputated Green 
function. 
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X 



[(7^7.7p^Lr") ® {YYi'PlT'') - 16(7^PlT") ® (t^PlT^)]} + O(e'). (70) 



c 



eff 



The Dirac structure in the last hne of the above equation vanishes in four dimensions. 
However, there is no way to express it as ex (simpler structure). The coefficient at this 
structure will give us the Wilson coefficient of an evanescent operator in the effective theory 
p4| . The necessity of recovering this coefficient (as well as keeping 0{e) parts of other one- 
loop coefficients) is a price we have to pay for regularizing infrared divergences dimensionally. 

The above result is the last one we need to know on the SM side. In the next subsection, 
we shall study the same Green functions in the effective theory framework. 

5.3. The effective theory side 

The lagrangian of the effective theory has been given in eq. (|l]). At present, we need to 
include in addition several non-physical operators. We write 

AC I 

C^cn.QEn{u,d,S,C,b,e,fl,T) + —l{ ^ V*,Vq,{C^,P^ + C^,P^ + C^.P.i) 

V^ [q=u,c 

+ E[(K*.K. + V:,V,,)Q + V,lV,,Cl]p] . (71) 

The operators P^ and Pi entering the effective lagrangian can be divided into three 
classes: physical, evanescent (i.e. algebraically vanishing in four dimensions) and EOM- 
vanishing (i.e. vanishing by the QCDxQED equations of motion, up to a total derivative). 

The physical operators have already been given in eq. (Q). However, for the purpose of 
the present section, it is convenient to redefine Pg so that it contains a sum over all the light 
charged fermions / weighted by their electric charges Qf 

P, = --,{sLi,bL)Y.QAh'f)- (72) 

^ / 

Such a redefinition of Pg does not alter its Wilson coefficient at leading order in electroweak 
interactions. 

As far as the evanescent operators are concerned, only Pf^ from the appendix will be 
needed in the present section. 

The gauge-invariant EOM-vanishing operators can be chosen as 

^31 = -isLrT%L)D''G;, + P4, 
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^32 = -^mf,SL00bR, 

r 

P33 = ^SL000bL, 
9^ 



P34 


i 
9 . 


SL a^'T^LG^ - G%SLT'^a^''0bL 


+ Ps, 


P35 


ie 


Sl B a^^'bLF^, - F^,7sLa^''0bL 


+ P7, 


-P36 


= -,{-SLrbL)d'F,, - Pg. 


Our sign convention in 


the covariant derivative acting on a quark field ip is 



(73) 



D,^ = {d^ + i9GlT^ + leQ^A^) ^. (74) 

The EOM-vanishing operators in eq. ( |73D can be assumed to contain the background 
gluon field only, because nothing but their tree- level matrix elements will be needed for the 
off-shell matching in the next subsection. However, a systematic off-shell renormalization of 
the effective theory requires introducing EOM-vanishing operators that contain the quantum 
gluon field as well. The explicit form of such operators is irrelevant here. Nevertheless, one 
should not forget that all of them enter into the sums over operators, such as the one in the 
last term of eq. (|7TD . 

It is not completely trivial to convince oneself that eq. (|73|) indeed contains all the gauge- 
invariant EOM-vanishing operators that we may encounter. One way to do this is to first 
write all the AB = —AS = 1 operators of dimension 5 and 6 containing the left-handed 
s-quark field only.0 The derivatives acting on the s-quark field can be removed by parts. 
One can start from writing down the 6 possible operators that contain the chromomagnetic 
and electromagnetic field strength tensors or their duals 

isLT^a^'bR)G%, isLT^rbL)D^G%, {-SLT-YD^bL)Gl,, 

Nothing new is obtained from the first two pairs of operators above, when the field strength 
tensors are replaced by their duals, because of the Bianchi identity and aajsl^ ~ Sap-^s^^'^^ ■ 
On the other hand, replacing the dual tensors by ordinary ones in the last pair of operators 
would break CP combined with 6 <-* s interchange even for mi, = and real CKM angles. 



^ Here, the dimension of an operator is understood as the sum of dimensions of the fields and derivatives 
it contains. Exphcit mass factors in the normahzation are not counted. 
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The remaining operators (apart from the four-fermion ones) must contain covariant 
derivatives. Since commutators of the covariant derivatives give field strength tensors, only 
one additional operator with three covariant derivatives (e.g. sl^D'^Bl) and one operator 
with two covariant derivatives (e.g. slD%l) remains. At this point, one has at hand a com- 
plete set of 8 gauge-invariant operators (apart from the four-fermion ones). The "magnetic 
moment" operators P7, P^ and the EOM-vanishing operators P31, ..., P36 are just certain 
linear combinations of them, P4 and Pg (up to total derivatives). 

Since both the u- and c-quarks are treated as massless in the present calculation, the 
lagrangian is symmetric under m <->■ c exchange. This symmetry has already been taken 
into account in eq. (|7TD : the same Wilson coefficients C^ occur both in the w-quark and the 
c-quark sectors. 

The lagrangian (|7lD is written in terms of bare fields and parameters. In order to express 
it in terms of the QCD-renormalized quantities, we replace 



ZgQ, TJib^Zmmb, ip ^ Zl^^ip, C^ ^^C^Z^i, (76) 



3 
for the QCD gauge coupling, 6-quark mass, quark fields and the Wilson coefficients, respec- 
tively. As far as the background gluon field G^^ is concerned, we only need to remember 
that gG^^ does not get renormalized. 

After QCD renormalization, the structure of the effective lagrangian is the same as in 
eq. (|7ID, but the Wilson coefficients G^ are replaced by some other constants that we denote 
here by A^ . Below, we shall need 



A'j = ZlJ2C?Zi^ forj = 1,2,4,11, 

Zm 2^ Q Zii + [Zm — 1) 2^ Q ^i(35) 
i i . 

Zm 2^ Q Zi^ + [Zm — 1) 2^ Cj ^j(34) 
i i 

^^ = Z^Z-'Y.C?Z,,. (77) 



^^ = z^zf 



A'i = z^,z;' 



For simplicity, we shall use the MS scheme in the present section. The MS results for the 
Wilson coefficients will be obtained later from the MS ones by simply setting 7^; — ln(47r) 
to zero, i.e. replacing k by ln(M^//ig). 
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In the MS scheme, the renormahzation constants read 
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for 5 active flavours, 
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The finite terms a^ can be different from zero if and only if P^ is an evanescent operator 
and Pj is not. Values of a?^ are fixed by requiring that renormalized matrix elements of 



evanescent operators vanish in 4 dimensions ||2J]. This requirement is just an extension of 
the MS'-scheme definition to situations where evanescent operators are present. 

Our off-shell operator basis is chosen in such a manner that as many operators as possible 
are EOM-vanishing. This means that no linear combination of the remaining operators is 
EOM-vanishing. In such a case, the EOM-vanishing operators do not mix into the remaining 
ones, i.e. Z^j = when Pj is EOM-vanishing and Pj is not. In consequence, we shall need 
to know explicitly only the mixing among the physical and evanescent operators.|^ 

The powers of coupling constants in front of our operators have been chosen in such a 
way that terms of order g'^'"' in the renormahzation constants originate from n-loop diagrams 
in the effective theory. As one can see, the sum of powers of gauge coupling constants in 
front of a given operator is always equal to "(number of fields in this operator)-4" . In the 
original QCD and QED lagrangians, the powers of coupling constants are equal to "(number 
of fields)-2" . Here, two powers are traded for Gp that normalizes the effective lagrangian. 

The renormahzation constants are found by calculating ultraviolet divergent parts of 
Feynman diagrams in the effective theory. When doing this, it is essential to clearly separate 
ultraviolet and infrared divergences. In order to do so, one can introduce an auxiliary mass 
parameter into all the propagator denominators (including the gluon ones), as explained in 
ref. [^. All the renormahzation constants in the effective theory up to two loops are known 
from the former anomalous dimension computations 0, §, ^, Q (although some of them 
need to be transformed to the "new" operator basis (Ef)). Here, we shall need the one-loop 
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Getting rid of ^^(34) and ^^(35), which enter eq. (|77|), is somewhat tricky - see subsection 5.4. 
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ai 



renorinalization constant matrix d^^ for {Pi, P2, P4, P7, Ps, Pg, Pu} only. It reads 

* 

1 



, (79) 





* 

where stars denote non- vanishing entries that are irrelevant for us. 

In addition, for the two-loop matching of photonic penguins in the charm sector, we shall 
need 

12 _ 116 ^22 _ n „01 _ n 
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16 
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16 
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,12 



116 

81 ' 

776 
243' 



,22 



^29 ~ 243' ^29 ~ 8l" ' '^(11)9 



148 
81 ' 



,01 



64 
27' 



At this point, we are ready to calculate all the necessary IPI Green functions on the 
effective theory side. This turns out to be very simple, because all the particles in the ef- 
fective theory are massless in our approach.]^ Consequently, all the loop diagrams vanish 
in dimensional regularization, because of the cancellation between ultraviolet and infrared 
divergences. In effect, we need to know only the tree- level matrix element of the effective 
lagrangian. The ultraviolet counterterms present in this matrix element reproduce precisely 
the infrared divergences in the effective theory, which have to be equal to the infrared diver- 
gences on the SM side. As we shall see, all the 1/e" poles will indeed cancel in the matching 
condition. 

External gluons in the Green functions considered on the Standard Model side have been 
the background ones. Therefore, we can maintain only the background gluon field in Ceff, 
since only tree-level diagrams are non-vanishing on the effective theory side. This is why we 
could omit EOM-vanishing operators proportional to quantum gluons in our operator basis, 
even though the calculation is performed off-shell. 

We now write down the effective theory counterparts of the Green functions considered 
in subsection 5.2. Their structure follows directly from tree-level Feynman rules for the 
operators given in eqs. (|^) and (|73|). 



^^ Remember that the 6-quark mass is formally treated here as a perturbative interaction with an external 
scalar field, and we include only terms that are linear in this interaction. 
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The b ^ S'j function reads (cf. eq. (|6lD) 

AC pP { ^^ ^^ 

V2 g [ j=i j=i 

with the coefficients at the structures 5*2, 5*8 and 6*10 given by 

h2 = -4A§5, hs = 2A^5-A^6, ho = A^, + Al,, 

h = -4A*5, h = 2Ai,-Ai,, /lo = A\ + Ai,. 

to all orders in QCD. Similarly, for b ^ s gluon we get 

12 12 



^2) 



AGfPrT'^ 
V2 ^ 



(k;k. + k:k6)E^a- + ytWtbT.'^^s,) (83) 



with 



M2 = 


= -4A^4, 


Ms = 


- 2A34 ^3]^, 


MlO = 


- A^ + ^34, 


W2 = 


= -4A^4, 


V8 = 


- 94* _ 4* 


^^10 = 


= Ag + ^134. 



(84) 

In both the b —>■ s'~f and b ^ s gluon cases, the coefficients at other structures depend on 
^32 and y433, too. In each of these two cases, coefficients at 12 independent Dirac structures 
Sj are given by linear combinations of only 6 independent quantities. It is just a consequence 
of QCDxQED gauge invariance of our effective lagrangian. Therefore, the coefficients at the 
structures S^ must satisfy 12 — 6 = 6 linear constraints. This must be the case also for the 
SM Green functions, because they must match the effective theory ones. Checking these 
constraints on the SM side has been an important cross-check in our calculation. 

The last function we have to consider on the effective theory side is the b — > sec one. It 
takes the form 

^^V:,V,,{A'i{^,PLn^{YPLn + A^,{^,Pl)^{YPl) 
+ Al, [{i.i.IpPlT'^) ® {^YYPlT'') - 16(7^PlT'^) ® {i^PlT"")]} 
+ [terms proportional to (^3^^ + A^)]. (85) 

5.4. The matching 

The Wilson coefficients can be perturbatively expanded as in eq. (|]). We shall ffist recover 
the Wilson coefficients at all the EOM-non-vanishing operators up to one loop. Then, two- 
loop contributions to the coefficients at P-j and Pg will be found. 
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A careful reader might be surprised that we start the matching without having considered 
diagrams with UV counterterms on the SM side. Apart from the electroweak counterterm 
proportional to s0b, we should include QCD renormalization of the quark wave functions 
and masses. 

The electroweak counterterm proportional to sPb is taken in the MOM scheme, at q^ = 
for the s0b term, and at vanishing external momenta for the terms containing gauge 
bosons. It is achieved by an appropriate flavour-off- diagonal renormalization of the quark 
wave functions. The only effect of such a renormalization in the present case is that the 
coefficients at the structure S'13 in eqs. (0), ( |B^ ) and ( p^ ) are completely renormalized away. 
This is welcome, because the structure S'13 was absent from the effective theory counterparts 
of these equations (eqs. (pl|) and (]83D). 

As far as the QCD renormalization of the quark wave functions in internal lines and in 
vertices is concerned, it combines to an overall factor, which could be obtained by renor- 
malizing only those terms in the vertices that correspond to external fields in a given Green 
function. However, one-loop external quark field renormalization is the same on the full and 
effective theory sides. Consequently, we can omit counterterms with Z^ on the SM side and 
simultaneously set Z^ to unity on the effective theory side. 

The same refers to the renormalization of the 6-quark mass, since rrib is actually treated 
as an external scalar field. We omit the corresponding counterterms on the full theory side 
and simultaneously set Zm to unity on the effective theory side. This is how we get rid of 
terms proportional to {Zm — 1) in eq. (ffTl). 

As far as the renormalization of the QCD gauge coupling is concerned, no such counter- 
terms occur on the full theory side in our particular calculation. On the effective theory side, 
we maintain all the necessary factors of Zg. 

The last relevant quantity that acquires QCD renormalization on the full theory side is 
the top quark mass. However, contributions from the corresponding counterterm diagrams 
can be obtained by differentiating lower order results with respect to nit (see below). 

Let us first match the b —>■ sec Green function up to one loop. The first thing to notice 
is that terms proportional to A^i + A^ in the last line of eq. (P3| ) are not important at the 
considered order, because 

A'^ = -A^, + Oig^). (86) 
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The reason for this relation is that the b -^ sdd IPI Green function acquires its leading 
contribution only at two loops in the SM. Lower-order tree-level contributions to this function 
must vanish in the effective theory, which implies the above relation. 

Similarly, from the fact that the h -^ se^e~ IPI function vanishes at one loop, we find 

A^ = +Ag + 0{g'), (87) 

so long as the VT-boson boxes and Z-boson penguins are not taken into account on the SM 
side (as we have assumed at the very beginning of this section). 

Returning to the b — > sec function, we compare eqs. (p9|), ( [70D and (0), and immediately 
find 



Al = T^AtW ---15-^0 +Oig\e\ 



^iV« (-^- - 15 - ^' 
(47r)2^^^ [ e ^^ 2 

A^ = -1 + 0{g' 



n 



^n = (1^(1-^-) (-7-^) +0{9\e), (88) 

which implies that (cf. eqs. ([771)- ([T^) with Z^ set to unity) 

Ci^(°) = 0, Cf) = -1, Cif)=0, (89) 

and 



^c(l) _ 






= 


/ 39 1 \ 
= -15 + 6fi; + e -— + 15k - 3k^ " 7^^^ ) + C'(e^), 

\ Zi Zi / 


(90) 


cf ) = 


= 0, 


(91) 


Oil 


= (1 - -) {-I -I)- ^Cf^all,, + 0{e) 




= 


- -\ + ^ + ^(^)- 


(92) 



Indeed, all the 1/e poles have cancelled in the final results for the one- loop Wilson coefficients. 

The coefficient C| is the only one that acquires a tree-level contribution in our calculation. 
For all the other coefficients considered below, we have Cj = 0. 

Let us now turn to the b ^ s gluon matching. Comparing eqs. (|66|)p| and (^31), and 



^^ Without 5i3, since it has been renormalized away by the electroweak counterterm mentioned in the 
beginning of this subsection. 
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solving the trivial set of linear equations |(p^,(]86|)|, one finds 



A'i 



^8 



^ivi" (i4" + 4" 



4 2\ 



(47r 



0(g\f- 



10 



0{g\e% 



(93) 



which implies that (cf. eqs. ([771)- ([7^)) 



cf) 



9 3 V54 9 3 18 / ^ ^ 



c: 



■c(l) 



l^\u 



(1) , „,(!) 
2 

11 1 



'l-eK)[-u\'> + ur^]+0{t' 



K + 0(e^ 



(94) 



Similarly, 



d« 



'1 — e/t 



a 



t(i) 



)Qt;«(x)+.;ff(a;))+0(a. 



(95) 



Finally, we perform the b —>■ s'j matching. Comparing eqs. (|6TD, (64) and (SI), and solving 
the trivial set of linear equations |(^),(^7|)|, one finds 



Ar 



9' 



{Any 

,4 



+ 



9 



(47r)4 

,2 



;i - 2en) {h!i^ + MJ) + 0{e\ 



+ 0{g' 



Al 



9 



+ 



(47r)2 
9' 



-h 



(1) 



h 



(1) 



2 --8 J+Oie' 



(47r)4 
which implies that (cf. eqs. ( [77| 



C 



c(l) 






iV(i) 

ATf (-!/.?) -/.?)) +0(e)]+0(,^), 
-(^) with Z^ and Z^ set to unity) 

(l-eK)(l/.« + /.S)+0(e^) 

23 /145 23 \ ,^, 2, 
36+4216-36'^)+^^^)' 



(96) 



38 4 

K + e 

27 9 



247 38 2 2 1 2 \ ^^/ 2 
162 27 9 27 ' ^ 



(97) 
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and 



C, 



42) 



;i - 2eK) Q/if + h 



[1 - 2e/t) 



112 107 



^e 243 
16 /23 145e 23e 



1 

e 

116 
8k 



I) 



3e V36 216 
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16 /I lie 



9e V3 1^ 



/€ +0(e) 



713 4 ^^, , 

— /t + C(e), 



243 81 



r42) 



NP{-lh?^-h^^^)-^^{all + Poall)C', 



u\ r^c(o) 



1 



12^c(0) ^llr~''=(l) I ^11/^^(1) 



'^29^2 ~l~ '^19^1 ~l~ '^49^< 



^01 ^c(l) 
"(11)9'-"11 



0(e) 



1 - 2e/t + e' 
128 



81e2 
16 



-1 



776 



+ 2/t' 



128 1496 5924 128 



^le2 243e 729 243 



TT 



— 

243e ' ' '' ^ 27e 



39 

-15 + 6/t + e ( h 15/t - 3/t^ TT^ 



27e 



7 2 /77 7 1 2 1 2 

1 K, + e K K, TT 

9 3 V54 9 3 18 



64/3 , ^^, , 



524 16 128 2 128 2 .^^ x 
729 3 81 243 ^ ^ 

Again, all the 1/e" poles have cancelled in the final results. 

Similarly, in the top sector we find 



t{i) 



c; 



'^9 



t{2) 



c;' 



(jtm 



[I — eK, 
[1 — en 
[1 - 2e/t) 
'1 - 2e/t) 



■(1) 



W, 



fri^) + no'i^) 



-lA'\^)-n"i^) 



cW, 






;ff\x) + ft\x) 



d 






1 r 



(9^ 



(99) 



^t(i) 



{a)\ + Po)Cr + ^7^8^ + 0{e) 



ii-^t(i) 



.(2) 



fr{x)-fr\x) 



(2), 



--I'l^^i-Cl''' --a'^Cf^ ^0{e). (100) 



e ■■" 9a; ^ e 

Here, the x-derivative terms stand for contributions from the top-quark mass renormaliza- 
tion on the full theory side. Instead of including these terms, we could just calculate the 
corresponding one-loop SM diagrams with counterterm insertions. However, derivatives give 
us the same results much faster. 

It is easy to verify that all the 1/e poles indeed cancel in C-f and Cg . As usual, 
the 0{t) parts of the one-loop Wilson coefficients have affected the results of the two-loop 
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matching. 

The results for Cf\ Cf\ Cf\ Cf ^\ C?^'\ C^^'\ C?^'^ and C^^'^ obtained in the 



present section have already been summarized in section ^, after passing to the MS scheme, 
i.e. replacing k, by ln(M^//iQ). All the other matching conditions summarized there have 
been found in an analogous manner. In the two-loop Z-penguin contributions to Cg and 
Cfo, the effect of renormalizing the spb term on the SM side was less trivial than in this 
section. In the two-loop matching for P^ and P2, some care was required at renormalizing 
the top-quark loop contributions in the MOM scheme. In addition, scalar integrals with 
three non- vanishing masses were necessary p3[. Nevertheless, the basic algorithm remained 



the same as in the Pj and Pg cases, which we have described in detail here. 

Summary 

We have evaluated two-loop matching conditions for all the operators relevant to 
B -^ Xsl^l~ in the SM. Details of this calculation have been presented only for the operator 
P7 and for the photonic penguin contribution to the operator Pg. As far as the remaining 
matching conditions are concerned, only the final results have been given. However, the 
method of the calculation was very similar in all the considered cases. 

Our results allowed to remove an important (~ ±16%) uncertainty due to the matching 
scale /io from the prediction for BR[B -^ Xgl'^l^] for low invariant mass of the emitted 
lepton pair {s G [0.05,0.25]). The obtained Standard Model prediction for the branching 
ratio integrated over this domain is 1.46 x 10^^. This result would change to 2.92 x 10^^ if the 
Wilson coefficient Cj (fj^b) had an opposite sign, as it might happen in certain extensions 
of the SM. 

There remains a sizeable (~ ±13%) perturbative uncertainty in the above SM result, 
which is due to the unknown two-loop matrix elements of the four-quark operators. Calcu- 
lable non-perturbative effects which have been included in our result are smaller than this 
uncertainty. Estimates of other non-perturbative effects suggest that they are not larger. 
Therefore, the next step in improving the accuracy of the theoretical prediction should be 
a calculation of the two-loop matrix elements of the four-fermion operators and one-loop 
matrix elements of the "magnetic moment" ones. 
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Appendix 

Here, we give the eight evanescent operators that were used in evaluating the anomalous 
dimension matrices given in section 0. Their explicit form defines what the MS scheme 
means in the effective theory. As before, the symbol Q stands either for u or for c. 

Pi, = (sL7Mi7M2 7M37M47M5&L)E(^7^'7'^^7'^'7^*7'^^g)-20P5 + 64P3, 

q 

Pi = {sLl,,^,,^,,^,,j,,T^QL){QLY'Y'Y'Y'Y'T%L)-20Pfi-256Pf^, (101) 

P22 = {-SLl,,l,,l,,l,,l,,TQL){QLr'r'r'l''l''hL)-2QPg-2hQP?, 

P25 = (sl7mi7m2 7m37m47m5 7m6 7m7&l) Y^^q^^^^^^^^^^^'^^^^^^^^^^^'^^q) - 336P, + I28OP3, 

q 

P26 = {sLl,,l,,l,,l,,l,,l^,l^rT''bL) J2{qY'Y'Y^^Y'Y'Y'Y'T'^q) - 336P6 + 128OP4. 

g 
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